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Quark-gluon vertex with an off-shell 0(a)-improved chiral fermion action 
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We perform a study the quark-gluon vertex function with a quenched Wilson gauge action and 
a variety of fermion actions. These include the domain wall fermion action (with exponentially 
accurate chiral symmetry) and the Wilson clover action both with the non-perturbatively improved 
clover coefficient as well as with a number of different values for this coefficient. We find that the 
domain wall vertex function behaves very well in the large momentum transfer region. The off-shell 
vertex function for the on-shell improved clover class of actions does not behave as well as the 
domain wall case and, surprisingly, shows only a weak dependence on the clover coefficient csw for 
all components of its Dirac decomposition and across all momenta. Including off-shell improvement 
rotations for the clover fields can make this action yield results consistent with those from the 
domain wall approach, as well as helping to determine the off-shell improved coefficient c,. 
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I. INTRODUCTION 

The study of gauge-fixed quark and gluon correlation 
functions on the lattice is a useful tool to understand the 
non-perturbative behavior of QCD. Within lattice QCD 
one may hope that the quark-gluon vertex in particular 
may be a momentum dependent probe that displays dis- 
cretization effects beyond those seen in the free field case. 
Indeed one may hope that knowledge of the quark-gluon 
vertex allows the determination of improvement coeffi- 
cients without resort to Ward identities. This may allow 
the non-pcrturbative 0(a) improvement obtained in the 
light quark limit for clover fermions ij using the PCAC 
relation to be extended to the Fermilab action |3| at ar- 
bitrary quark mass. 

Beyond lattice QCD knowledge of the non- 
perturbative n-point functions gives the only theoretical 
access to the deeply infra-red sector of the theory and 
may be used to provide non-perturbative information 
for model calculations. Further we can obtain the 
non-perturbative running of a physically defined cou- 
pling, and may be able to use these non-perturbative 
correlation functions in the Dyson-Schwinger equations. 
DSEs are of qualitative and quantitative importance 
in understanding dynamical chiral symmetry breaking 
and confinement, and simulations of lattice QCD also 
provide access to QCD's Schwinger functions. Previous 
studies of the quwiched theory have yielded gluonQ 
and clover quark 4] propagators and the quark-gluon 
vertex on the lattice has also been studied extensivelyj5|- 
These previous studies of the vertex were limited to 
the on-shell improved, non-chirally symmetric clover 
action, and post-simulation field rotations were used to 
implement off-shell improvement. 
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The domain wall chiral fermion action [g is automat- 
ically 0(a)-improved, and domain wall fields should 
greatly aid studies of off-shell correlation functions. 
The off-shell quark propagator has been previously and 
comprehensively studied for the domain wall fermion 
action jj, and renormalization constants were determined 
for many fermion bilinears in the chiral limit for the Wil- 
son and DBW2 gauge actions. The off-shell properties of 
the propagator were found to be remarkably continuum- 
like for this action, and one expects that this would also 
be seen for the domain wall vertex function. 

In this paper we also use this more advanced domain 
wall fermion action, which is chirally symmetric (bar ex- 
ponentially small terms), and demonstrate that such im- 
provement is indeed seen in our calculations of the quark- 
quark-gluon three-point function and the vertex function. 

Further, one might expect that a Dirac structure de- 
composition of the vertex function could be used to 
match the clover action (with field rotations) to the do- 
main wall action. This could provide a mechanism for 
fixing the clover coefficient non-perturbatively without 
the imposition of the PCAC relation, and it was hoped 
that this would allow a non-perturbative clover coeffi- 
cient determination for all quark masses However we 
find weak dependence on this coefficient after absorbing 
a trivial mass renormalization and instead find the cor- 
relation function is most sensitive to the improvement 
parameter. 

The structure of this paper is as follows: Section lITIde- 
fines the notation and lattice correlation functions mea- 
sured to obtain gauge-fixed two-point and three-point 
fimctions. Section UTTI lists details of the simulation pa- 
rameters. Section IIVI presents the gauge-fixing strat- 
egy used. Sections |3 and IVII present simulation results 
of the quark and gluon two-point functions. Section 
IVIII presents simulation results of the quark-quark-gluon 
three-point function and the truncated vertex. Section 
IVIIII discusses the off-shell improvement and contrasts 
the DWF and clover cases. Finally, Section Hxl presents 
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a summary, conclusions, and future outlook for this pro- 
gram. 



II. NON-PERTURBATIVE GREEN 
FUNCTIONS: DEFINITION 

A. Off-shell 0{a) improvement 

The on-shell 0(a)-improved Wilson clover action|9| 
can be written as 
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To off-shell improve the physical quantities calculated 
from this Wilson clover action, we keep the action un- 
changed but perform post-simulation rotations'lO] to im- 
proved quark fields of the form 



ip — > tp 
ip — > Ip 



1 -I- aCq{-Dc +m) - acNGi7M ^i^ 



(4) 



The improved quark propagator Si{x) = {tp{x)'4'{Q)) 
can then be expressed in terms of the unimproved quark 
propagator. In momentum space, we have: 

Slip) = Zq [{SL{p) + 2ac'^) 

+ ac^Gi\ili-i^,SL{p)^ +0{a^) (5) 

where 5*^ is the unimproved quark propagator. The con- 
stant cngi has been previously demonstrated to be zero 
at tree levelj^ and small for light quark masses[lT|. 
The domain wall fermion (DWF) action is 
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(1 - 75)(5s+l,s' + (1 + l5)5s-l,s' - 2(5s,s 
(1 - 75)<5s,L,-1^0,s' + (1 + 75)<5s,0<5l,-1 



(8) 

where s and s' lie in the range < s, s' < — 1, Afs is 
the five-dimensional mass, and m/ directly couples the 
two domain walls at s = and s = Lg — 1. The DWF 
action is 0{a) off-shell improved due to its exponentially 
accurate chiral symmetry 'T^, and no further improvement 
in the action or quark fields is performed. 



B. The gluon field 

The lattice gluon field A"^ is defined by 
1 



Al{x + fi/2)T- = 



2igo 



{U^{x)-Ul{x)) 



^TT{U^ix)-Ulix)) 



0(a3),(9) 



(7) 



where a is the color index, /i is the Lorentz index, go is 
the bare coupling constant, and Trlr^r''} = ^Sab- 
The gluon propagator (in Landau gauge) is 

D;1{q) - :^(A^(<z)A^(-g)) ^ SabP,.D{q'), (10) 

where D{q^) - ^ E^,a DZiq) and P^, = S^. - 



C. The quark-gluon vertex function 

The quark-gluon vertex function can be calculated 
from the following three-point function^: 

V^^{x,y,z)% = {Si^^ix,z)A';,{y)). (11) 

where the lattice gluon field A° is defined in Eq. IHl 
The Fourier-transformed, amputated and transverse- 
projected vertex function in momentum space for Landau 
gauge is given by 

A^'''''^'{p,q)% = A,{p,qUT^, - P,.A:''^'{p,q)% 

^{S{p)r'P,,V^\p,q)%{S{p + q))-'{D{q))-\l2) 

where the projected vertex function is introduced since 
the gluon propagator is singular for general momenta, 
posing problems for a full amputation. Note that the S 
means Sl{p) for the bare vertex and Si{p) for the off- 
shell improved clover quark field. 

In the case of g = the gluon propagator is no longer 
proportional to the projection P^jy, but rather is 

Kf^\p,q^O)% 

^{S{p))-'P,,V,^{p,q^0)%{S{p))-'{D{0))-'{13) 

The non-perturbative gluon propagator in a finite volume 
remains well defined and the same amputation process is 
used. 
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III. SIMULATION PARAMETERS 

All of the calculations in this paper were performed 
with the Wilson gauge action^^l the quenched approx- 
imation at /3 = 6.0 , = 1.922(40) GeV with a 16^ x 32 
lattice size 01 . The simulation was performed using a 
heatbath algorithm with 2000 thermalization sweeps and 
500 measurements, separated by 2000 sweeps each. 

The fermion actions used are DWF and the Wilson 
clover action with various parameters. The details are 
listed in Table 0] and we note that the pseudoscalar 
masses for three of the clover action combinations were 
(somewhat approximately) matched to the domain wall 
pseudoscalar mass using a single, low-statistics pass. 



two ways. The first is naively fixing to Landau gauge 
in the usual way, and the second involves fixing to Lan- 
dau gauge from a deterministic starting gauge orbit, by 
first fixing to maximal axial gauge 'l 55 . This is a common 
approach which gives a deterministic Gribov copy, inde- 
pendent of the gauge of the initial lattice configuration. 

We found that the two different procedures agree for 
the gluon propagator, have slight effects on quark prop- 
agators (decreasing as the momentum increases) but 
give very different results in the vertex calculationjSj, as 
shown in Figure Q We believe the differences are caused 
by Gribov copies and from now on, we will follow the de- 
terministic procedure for better control over this effect. 



TABLE I: Simulation parameters 



Label 


Parameters 


PS meson 


DWF 


Ms = 1.8, Ls = 16, 
ma = 0.0125 


442(9) MeV 


csw ~ 1.769 
csw = 1.5 
csw = 2.0 

csw = 1.479 


K = 0.1346, csw = 1.769 
K = 0.1382, csw = 1.5 
K = 01313, Csw = 2.0 

K = 0.1370, csw = 1.479 


470(10) MeV 
462(11) MeV 
452(13) MeV 
796(8) MeV 



We use periodic boundary conditions for both the 
spatial and temporal directions to calculate both the 
fermionic and gluonic propagators. The momenta used 



for Fourier transformation are 

I e {0,1,2} and \nt\ G {1,2}. 



27rn^/L^, with 




(pa) 



FIG. 1: Comparison between Landau gauge fixing with 
(squares) and without (diamonds) first fixing to maximal ax- 
ial gauge for the vertex function. 



IV. GAUGE FIXING STRATEGY 

Since we are interested in gauge-dependent off-shell 
quantities we choose to fix to a convenient gauge, namely 
Landau gauge, by maximizing the following functional: 

F[U^{n),g{n)] = 1 - ^ ReTrC/9(n), (14) 

where Uf^{n) = g{n)Uf^{n)g{n + /t)^ with gauge transfor- 
mation matrix g{n). The gauge-fixing stopping condition 
is 

^^Tr(B(n)i?(n)t) <10-« (15) 

where N is the number of total sites on the lattice and 
B{n) = C{n) - C(n)t - ^Tr (C(n) - C(n)t) , (16) 
C{n) = (c/^«(n)t + C/^^(n-A^)). (17) 

Landau gauge suffers from Gribov ambiguities 0|, and 
due care must be taken to ensure our results are not sig- 
nificantly affected. We checked for the influence of Gri- 
bov copies on each of our observables by gauge fixing in 



V. FERMIONIC TWO-POINT GREEN 
FUNCTION 

We start from the two-point function used for non- 
perturbative renormalization in the RI-MOM scheme 1^. 
The quark and gluon propagators from quenched lattice 
configurations have been well studied previously for both 
fermion actions^H- However, these remain a prerequi- 
site for our exploratory study of the vertex function, and 
our results must be defined and introduced prior to mov- 
ing on to the vertex function. 

The bare quark field renormalization factor, Zq(p^), is 
calculated from 

^^(^")- 12E>.a)- ^^(^"-'""-^^")' ^''^ 

where, as we shall see, constraints on the above quantity 
in the matching window allow Zq to be determined from 
a fit. 

In Figure 13 we show the bare quark field renormaliza- 
tion Zq{{paf') for each of the actions listed in Table HI 
using the bare fields. We calculate the Zq factor in the 
RI-MOM scheme by removing the expected perturbative 
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FIG. 2: The quark field renormalization factor Zq for differ- 
ent actions. 



FIG. 3: Tlie quarlc field renormalization factor Zq for the 
bare, SI and MS schemes on our DWF data set. 



momentum dependence to define a scale invariant quan- 
tity 



(19) 



and then fit the (ap)^ dependence of this SI expression to 
the form /i + /2(pa)^ where the /2 term absorbs lattice 
artifacts. 

From this we obtain the phcnomcnologically relevant 
renormalization factor between the lattice and the MS 
scheme. 



^RI 



in 



^(l)RI^ 



at 



(4^) 



2^0 



(2)RI 



(20) 



where definition and values for CA((pa)^), ctg, Zq^^^^ and 

Zq^^^^ are well known 7] . Figure O shows Zq for these 
three schemes for domain wall fermions. 

Table Ull shows our results for the various fermion ac- 
tions under consideration. As one would expect, where 
previous data exists for the same action, these results are 
consistent. In particular, for the domain wall action we 
refer the reader to Ref. where the massless limit was 
taken. 



TABLE II: Quark field renormalization factor for the fitting 
range of 0.8 < {paf < 2.0. 



Label 


^si 


^MS 


DWF 


0.802(9) 


0.806(9) 


csw = 1.479 


0.825(6) 


0.829(6) 


csw = 1.5 


0.840(6) 


0.844(6) 


Csw ~ 1.769 


0.846(6) 


0.850(6) 


csw = 2.0 


0.849(6) 


0.854(6) 



Following Ref. [l^, we expect that in the large mo- 
mentum region the scalar part of the quark propagator 



o 


DWF 


□ 


c_^= 1.479 


o 


0^^=1.5 


A 


c_^=1.769 


X 





1.5 2 
P (OeV) 



FIG. 4: The quark field renormalization factor Z^^{{pa) 
for different actions. 



behaves like 



Y^Tr(5i(p)) 



C X (pa)2 + {^2c'g + 2cNGiZg) 

+^^ + o(l). 



(21) 



As mentioned in Subsection lll Al cngi is zero at tree level 
and has been found to be numerically very small in the 
light quark region'll], enabling us to take it as zero in 
our calculation. Therefore, we can use this formula to 
determine c^. 

The resulting parameters, used later in this paper, are 
displayed in Table UTTl We note that the first two terms 
in Eq. 1211 violate chrial symmetry and are not allowed 
in the DWF case. Our values for the clover actions are 
plausible given tree-level estimates of around 0.25. 
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TABLE III: Fitting to the first three terms in Eq. ^ with 
fitting range of 0.8 < (pa)^ < 2.0 produces the improvement 
coefficient c' . 



Label 


/ 


csw = 1.479 


-0.247(19) 


csw = 1.5 


-0.258(19) 


csw = 1.769 


-0.257(19) 


csw = 2.0 


-0.252(19) 



VI. 



GLUONIC TWO-POINT GREEN 
FUNCTION 



0.5 



1.5 

q (GeV) 



2.5 



There are many renormalization schemes that may be 
apphed to the gluon propagator, and one of the more 
common in the hteraturc is the MOM scheme where the 
gluon field renormalization factor is constrained by re- 
quiring that, at some fixed scale /i, an effective tree level 
Coulombic form holds: 



Z^\t,)D{q';t,)\ 



1 



(22) 



The quantity ^3(5^)^^ = q^D{q^) is displayed in Fig- 
ure[Sl and we see that this condition leads to a renormal- 
ization constant that is highly sensitive to the scale /i 
in the momentum regions currently accessible to lattice 
QCD. 

This scale dependence is an artifact of the renormaliza- 
tion condition, and we pick a fixed renormalization point 
a/i = 1 and use the one-loop form 17] of the running to 
present the data in a manner that should be perturba- 
tively scale- invariant: 



{Zf{q^))-'^q-D{q-) 



(23) 



where do — ^ for the quenched approximation in Lan- 
dau gauge, a/i = 1 by choice, and A = 0.35(5) is obtained 
from a fit to the propagator. We note that the error is 
obtained from the sensitivity of the result to the upper 
limit of this fit range in the region 0.5 < (ag)^ < 2.0. 
Our value for Aa is identical to that obtained in Ref. 
when fitted over the range used in that paper, but it 
shows great sensitivity to the range. 

We can see that this scale invariant presentation in 
Figure shows a nice plateau demonstrating perturba- 
tive behavior and gives confidence in our result at and 
around this renormalization point. We obtain the MOM 
scheme renormalization constant as {Z^^)~^ = 2.45(1) at 
the scale of /i = a'^ = 1.922(40) GeV. 

A more conventional approach would be to use the MS 
scheme at some scale. However, the conversion to this 
scheme is not known and we isolate this lack of knowl- 
edge, localizing it to a single continuum renormalization 
constant Zy,{ji) which converts from the RI-MOM scheme 
at a/i = 1 to the MS scheme at momentum scale JL. 



FIG. 5: The naive RI-MOM scheme renormalization constant 
Z^^ — q^D(q^) in Landau gauge. 
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FIG. 6: Scale invariant MOM scheme renormalization con- 
stant Zf^(g^)~^ at fixed renormalization point (a/i)^ — 1. 



VII. QUARK-QUARK-GLUON THREE-POINT 
FUNCTION 



We consider the three-point function consisting of a 
gluon with incoming momentum q and two quarks with 
incoming momentum p and — {p+q) as shown in FigureH 



Quark-Gluon Vertex 




-> A^'''(p,q) 4- 



-ip-Ki) 



FIG. 7: Quark-quark-gluon momentum space, three-point 
function. 
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Lorentz structure of vertex 



D. Tree- level lattice perturbation theory 



In the continuum the vertex function has been ex- 
pressed in terms of its general Lorentz structure in 
Ref. fsj. This decomposition is complex and requires a 
comprehensive set of lattice momenta to resolve the dif- 
ferent components. 

In the following subsections we will study kinematic 
regions that simplify the form factors appearing in the 
decomposition allowed by Lorentz symmetry in each case. 
As usual, appropriate spin and momentum projections 
can be used to determine the form factors from the bare 
lattice data. 



B. Hard recoil kinematics 

We first require both quark lines to carry momentum p 
giving a gluon momentum of —2p. In this case the vertex 
function may be decomposed as follows 



Here we note that the tree-level vertex function for the 
clover lattice action (without off-shell improvements) is 



(24) 



The form factors \'i and r5 can be easily separated by 
different Lorentz projections. 



C. Soft gluon kinematics 

Next we choose the point g = 0. At this kinematic 
point the tensor contribution vanishes, and we may write 
these form factors in terms of their popular Lorentz de- 
composition as follows: 

- 2iX^[p')p^]. (25) 

The form factors Ai and A2 are coupled by the same 
Lorentz projections but can be decoupled as following: 



Ai(p') 

Hp') 



1 

48 (pa) 
i 



^)Alf*)(26) 

4, 



5:Tr(7X")-iAi(P^) 
24(,«).ETrM-). 



(27) 
(28) 



Our method projects Ai and A2 for general kinematics 
while in Ref. Q the authors use only kinematic points 
with = to isolate Ai from A2. 



a!"1' 



lat 



I sm I p^a + -q^a 



icsw cos ( ig^a j ^ a^,^ sin (q^a) . (29) 



Thus, we might expect that to leading-order the clover 
action will be most sensitive to the clover coefficient 
through terms proportional to cr^^. Therefore, in investi- 
gating the feasibility of fixing the clover coefficient from 
the vertex function we will consider the vertex function 
defined from bare (un-rotated) quark fields. We shall 
see later that, somewhat interestingly, this turned up a 
negative result. 

Note that in subsection IVII El we will use the "bare" 
clover action to calculate the various form factors to avoid 
contamination from the poor choice of off-shell improved 
coefficients, while in subsection IVII F1 we will focus on 
the off-shell improved vertex with the non-perturbatively 
improved clover action. 



E. Form factors 

Figures I8I12I display the lattice results for Ai(p^), 
A2(p^), Xsip^), Ai(p^), and r^^p"^) in the MS scheme us- 
ing the renormalization constants calculated earlier in 
this paper, up to the unknown constant Z^{jl) in Sec- 
tion This is shown for both the bare clover actions 
and the domain wall action. We will discuss the off-shell 
improvement of the clover action via post-simulation ro- 
tations in the following section. 

We note that the off-shell improved DWF action in 
general shows a marked difference in the A3 scalar form- 
factor; see Figure ^2 in particular. The difference be- 
tween the DWF and clover action results is likely to be 
caused by the off-shell improvement present with domain 
wall fermions, and one might hope that the form factor 
A3 would be a good means to adjust the off-shell improve- 
ment coefficients for clover actions. There will be more 
discussion in Section rVllI*'! on this possibility. 

The quantity X'^ is the form factor reflecting the 7^ 
term in the fermion action. This quantity appears to be 
the most similar among the actions studied here. 

As seen in Figure 1121 we find no statistically signifi- 
cant evidence of csw dependence in the a projection of 
the vertex function, despite both its leading order con- 
tribution to the tree-level vertex and our examining it at 
a sizeable range of lattice momenta. 

The difference between the (bare) clover and DWF ac- 
tions is again quite marked and may be in principle useful 
for determining off-shell improvement rotations. How- 
ever, since csw (in perturbation theory) couples with this 
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FIG. 8: The vector vertex function form factor Xi{p^)/ Z^lJt) 
for the soft gluon case. This is displayed for both the chirally 
symmetric DWF action and the bare clover action. 
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FIG. 9: The vector vertex function form factor X2{p'^) / Z^lJt) 
for the soft gluon case. This is displayed for both the chirally 
symmetric DWF action and the bare clover action. 



vertex function at leading order and since the non-zero 
gluon momenta are statistically noisier, we will only cal- 
culate the from A3. 

F. Determining c'q directly from the vertex function 

The quantity c'^ appearing in the improved quark field 
can be calculated from the quark pro pagator alone, as 
demonstrated previously in Refs. [lO, [llJ| and in this 
work. Here we propose and will demonstrate that we 
can determine this coefficient independently by match- 
ing form factors to the already improved DWF vertex 
function for the q = kinematics. Any consistency be- 
tween these two independent determinations of will 
contribute to the body of evidence that cngi is near zero. 



1. Tree-level calculation 

If we apply the off-shell quark field improvement of 
Eq.Q), we obtain the following expression for the tree- 
level vertex function 

^ {ism{pf,a + qf,a/2){~l + 2cg) 
+ icNGi [sin (p^a -|- g^a) -|- sin (p^^a)] cos (p^a -|- 
- {cos {p^a + q^a/2) {1 + AcqTua)} 
+ io'i^Li' {i/2 cos sin (q^a) (csw - 2c,) 

+ «CNGI [sin (p^a + q^a) - sin {puo)] cos (p^a + q^a/2)} 

(30) 

which depends on the coefficients of (= — 0.5cq) and 
Cngi- As above we assume cngi is zero, so the only task 
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FIG. 10: The scalar vertex function form factor X3{p^) / Z^lJt) 
for the soft gluon case. This is displayed for both the chirally 
symmetric DWF action and the bare clover action. 
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FIG. 11: The vector vertex function form factor A'j (p^) / Z:i{jt) 
for the hard-recoil kinematics. This is displayed for both the 
chirally symmetric DWF action and the bare clover action. 



is to determine the quantity from the vertex function. 

This expression illuminates how the different form fac- 
tors might be expected to depend (to leading order) on 
the different improvement coefficients. In particular, we 
note the leading-order sensitivity of the cr^,y piece, rs, to 
the clover coefficient. 



2. Numerical work 

Here we consider the clover action with the non- 
perturbative csw coefficient. We included the ma depen- 
dence (from on-shell improvement) by replacing the Zq 
in Eq.l^l with (1 + bqUiajZq for the clover action and 
performed the same calculation steps as in the unim- 
proved field cases. Since the ma is small and bq is 1.0 



at tree level, we expected the contribution from bq to be 
negligible. Indeed, we saw no significant changes in the 
calculation on quark propagator and form factors from 
quark-gluon vertex function. 

In Figure El we compare the three form factors Ai, 
A2 and A3 between the domain wall and bare clover ac- 
tions where the latter are computed for various values 
of c'q e {0.0, -0.20, -0.257} and with bq taking its tree- 
level value, bq = 1.0. As we predict from our tree-level 
exercises, the form factor A3 is indeed very sensitive to 
c'q. The calculation shows that our value for c'q calcu- 
lated from quark propagator is closest to the DWF re- 
sult, although at large momentum, the (pa)^ effects domi- 
nate and disagree with the completely off-shell improved 
result. Clearly we successfully confirm an independent 
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FIG. 12: The tensor vertex function form factor rs (p^) j Z-i (Jt) 
for the hard-recoil kinematics. This is displayed for both the 
chirally symmetric DWF action and the bare clover action. 



method to determine c^. 



VIII. DISCUSSION 

We have calculated the form factors Ai, A2, A3, A'j^ 
and Ts for both the bare clover action and the DWF 
action in two kinematic regions. We have numerically 
demonstrated rather clearly the automatic off-shell im- 
provement of chiral fermion formulations. 

Performing post-simulation rotations on the clover 
data can yield improved form factors. This involves wor- 
risome systematic errors, such as the taking of cngi — 0, 
which are entirely avoided in the DWF formulation. 

We demonstrate here that the use of the vertex func- 
tion form factors as a means of fixing at least one im- 




FIG. 13: The form factors Ai, A2, A3 (from top to bottom) 
calculated with q — kinematics and with various c'g. The 
various actions are represented by different symbols: DWF 
(circles); NP clover with c'g = 0.0 (squares), c'g = —0.20 (tri- 
angles) and c'g — —0.257 (crosses). 



provement coefficient of the clover action in a manner 
that does not require the light quark limit. Unfortu- 
nately we have not been able to fix the crucial on-shell 
coefficient cgw- 

From the vertex form factors for the q — kinematics, 
we showed that we can provide an independent way of 
determining the coefficient appearing in the improved 
quark field for the non-perturbatively improved clover 

= —0.257 is very close to the 
0.25 with the same tree-level 
-0.285 with bq = 1.14, 



action. Our finding of 
tree-level prediction Jilj, - 
bq, or to the mean-field value of 

and also to previous calculations'! 11 from the quark prop- 
agator alone. 

We note that simulations with different clover coeffi- 
cients and matched pseudoscalar masses produced iden- 
tical vertex function form factors. One would naively 
expect to see a difference due to different choices of 
clover coefficient, particularly in the T5 form factor and 



10 



in the large momentum region where the leading pertur- 
bative behavior should dominate and the universal cut- 
off-insensitive infrared physics disappears. However, it 
seems that the we have not yet reached a sufficiently fine 
lattice spacing where this perturbative intuition applies. 
Thus this coefficient cannot be fixed at these lattice spac- 
ings using this approach. 

In contrast to the insensitivity to the clover coefficient, 
we sec a large change in ofF-shcll behavior between the 
entire clover class of actions and the domain wall action. 
In particular, since the domain wall action is 0(a) off- 
shell improved and this (un-rotated) clover action is not, 
large differences serve to demonstrate the effectiveness of 
the domain wall action, and it surely is a better action 
to study the vertex function in the future. 



fixed pseudoscalar mass. Studying various Dirac struc- 
tures at large momentum for the domain wall action 
showed vast improvements in the off-shell behavior of the 
theory. Also promising is the use of the behavior of the 
vertex function as a new method of fixing the 0(a) off- 
shell improved field coefficients for otlic;! fcrmion actions, 
in particular during the non-perturbative step scaling of 
relativistic heavy quark actions. 

Comparison with the less-improved clover action al- 
lowed us to provide a new method to determine the coeffi- 
cient c'g through the form factor A3. The result found was 
consistent with the value calculated from the quark prop- 
agator and previous published results. Further, this con- 
sistency may be taken as additional evidence that cngi 
is negligible. 



IX. SUMMARY AND OUTLOOK 

We have studied the non-perturbative structure of 
two- and three-point functions for both the domain wall 
fermion action with an automatically 0(a) off-shell im- 
proved quark field, and for the on-shcU improved clover 
action with various choices of the csw coefficient at fixed 
pseudoscalar mass. For the clover action the vertex func- 
tion shows weak dependence on various choices of csw at 
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